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Abstract 

To meet the increasing demand of wood biomass worldwide in the context of climate change, developing improved 
forest tree varieties for high productivity in water-limited conditions is becoming a major issue. This involves breeding 
for genotypes combining high growth and moderate water loss and thus high water-use efficiency (WUE). The present 
work provides original data about the genetics of intrinsic WUE (the ratio between net C0 2 assimilation rate and sto- 
matal conductance, also estimated by carbon isotope composition of plant material; 6 13 C) and its relation to growth 
in Pinus pinaster Ait. First, heritability for 6 13 C was estimated (0.29) using a 15-year-old progeny trial (Landes prov- 
enance), with no significant differences among three sites contrasting in water availability. High intersite correlations 
(0.63-0.91) and significant but low genotype-environment interactions were detected. Secondly, the genetic architec- 
tures of 6 13 C and growth were studied in a three-generation inbred pedigree, introducing the genetic background of 
a more-drought-adapted parent (Corsican provenance), at ages of 2 years (greenhouse) and 9 years (plantation). One 
of the quantitative trait loci (QTLs) identified in the field experiment, explaining 67% of the phenotypic variance, was 
also found among the QTLs detected in the greenhouse experiment, where it colocalized with QTLs for intrinsic WUE 
and stomatal conductance. This work was able to show that higher WUE was not genetically linked to less growth, 
allowing thus genetic improvement of water use. As far as is known, the heritability and QTL effects estimated here 
are based on the highest number of genotypes measured to date. 

Keywords: Breeding, carbon isotope composition, genotypexenvironment interaction, genetic variability, growth, heritability, 
maritime pine, QTL, water-use efficiency. 



Introduction 

From 2000 to 2010, the world lost 130 million hectares of 
forest, corresponding to 3.2% of the total forest area, but 
on the other hand, it has gained 78 million hectares, mainly 



as planted forests and natural forest expansion (Food and 
Agriculture Organization, 2012). This increase in planted for- 
ests corresponds to an increase in wood demand, mainly for 
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fuel and industrial timbers, since natural forests cannot meet 
this demand anymore (Sutton, 1999) 

Climate change will affect planted forests, through 
increased tree mortality (Allen et al, 2010) and changes in 
net primary production (Hanson and Weltzin, 2000). While 
increasing temperature and rising atmospheric C0 2 concen- 
tration are likely to increase growth (Way and Oren, 2010; 
Penuelas et al, 2011), a reduction in precipitation and a 
change of its distribution in the annual cycle will result in 
increased soil water deficits during the growing period and 
therefore in a reduced biomass production. Based on a mod- 
elling approach, Loustau et al. (2005) predicted a decrease of 
net primary production for Pinus pinaster forests in south- 
western France because of a pronounced shift in seasonal 
rainfall from summer to winter. Dvorak (2012) argues that, 
considering the impact of plantation forestry on the hydro- 
logical cycle, the issue of water use will be one of the most 
important questions in many countries in the coming decades. 

Two complementary strategies are being developed to 
adapt planted forests to these new climatic conditions: adap- 
tive management of forest stands and adaptation of genetic 
material, either by planting new species which will use less 
water for biomass production and can better cope with a 
decrease of water availability or by adapting endogenous 
genetic resources by breeding. Existing forest tree breeding 
programmes mostly rely on improving wood biomass pro- 
duction and quality as well as disease and pest resistance 
(reviewed by Mullin et al, 2011 for conifers). Adaptation to 
climate and soil characteristics is generally considered indi- 
rectly, by evaluating growth or survival. Given the challenges 
imposed by climatic changes, it becomes necessary that forest 
tree breeding also includes selection criteria related to the use 
of water. 

To optimize the overall amount of water versus biomass 
accumulation, but also the maintenance of growth under 
reduced precipitation, is a crucial objective, which can be 
addressed by evaluating the ratio between biomass accumu- 
lation and water used (i.e. water-use efficiency, WUE). The 
direct estimation of WUE at the whole-tree level, for exam- 
ple by estimating time -integrated transpiration and biomass 
accumulation, is difficult (Kruse et al, 2012), especially 
for adult trees, and it is impracticable to phenotype large 
populations of plants for estimations of genetic param- 
eters. Screening of WUE of a large number of genotypes 
is commonly based on using a proxy, the carbon isotope 
composition (5 13 C) of plant organic material, such as nee- 
dles or wood (Brendel et al, 2002, 2008; Dillen et al, 2008; 
Gonzalez-Martinez et al, 2008). Indeed, there is a linear pos- 
itive relationship between carbon isotopic discrimination — a 
dimensionless measure of depletion in 13 C in plant material 
as compared with atmospheric C0 2 (carbon source for pho- 
tosynthesis) — and intrinsic water-use efficiency (W\) — the 
ratio of net C0 2 assimilation rate to stomatal conductance of 
water vapour (A/g s , Farquhar and Richards, 1984). For plant 
material grown in common environmental conditions, the 
S I3 C of atmospheric C0 2 is assumed common among plants 
and thus single measurements of 5 13 C of plant material can 
be used to estimate W v W{ represents the leaf-level processes 



involved in whole-plant WUE, and screening for differences 
in W { among tree genotypes is facilitated by the time- and/or 
crown-integrative nature of leaf or tree ring 6 13 C. 

The main goal of the presented research was to study the fea- 
sibility of integrating WUE as a new selection criterion in the 
maritime pine breeding programme, similarly to the success- 
ful approach applied in breeding wheat cultivars in Australia 
(Richards, 2006). Maritime pine (Pinus pinaster Ait.) was used 
for several reasons. First, a breeding programme was launched 
in France in the 1960s to improve growth and straightness, and 
genetic resources are available to work both on genetic param- 
eters and genetic architecture. Second, the planted forest in 
southwestern France, where this species spans over 1 million 
hectares, is expected to undergo reductions in water availabil- 
ity in the context of climate change, particularly during the 
growing season (Loustau et al, 2005); water availability has a 
major impact on maritime pine growth (Aranda et al, 2010; 
Corcuera et al, 2010; Sanchez-Salguero et al, 2012). Third, 
previous reports have shown that W { (as estimated from 6 I3 C) 
correlated well to WUE at the whole-plant level (Guehl et al, 
1995, 1996) and displayed a large natural variability (Guehl 
et al. , 1995; Correia et al. , 2008), which might partly be due to 
differences in W y among natural populations (Tognetti et al, 
2000). 6 13 C has also been shown to be heritable in this species, 
with low (Aranda et al, 2010) to medium heritability (Brendel 
et al. , 2002; Lamy et al, 201 1), although it has been generally 
estimated with a rather low sample size and depends on the 
experimental set up and the genetic background. The genetic 
architecture of WUE has already been studied in one mari- 
time pine family but, similarly to heritability, has also been 
estimated in a reduced genetic background and with a limited 
number of offspring (Brendel et al. , 2002). 

In this context, the specific objectives of this work were two- 
fold: (i) to estimate genetic parameters for 5 13 C and growth, 
based on the Landes maritime pine breeding population, using a 
progeny trial established in three sites with contrasting soil water 
availability; and (ii) to compare the genetic architecture of 6 13 C 
and growth across different tree ages in a cross between a geno- 
type from the breeding population and a genotype from a region 
with lower soil water availability. The rationale for the introduc- 
tion of 5 13 C as new selection criteria for maritime pine breeding 
is discussed with the aim of providing improved varieties that 
maintain biomass production under water-shortage conditions. 

Materials and methods 

Plant material 

Two types of plant material (Supplementary Table SI available at JXB 
online) were used to study the genetics of WUE and its relation to 
growth: (i) a progeny trial from the Landes breeding population to esti- 
mate genetic parameters; and (ii) two F2 quantitative trait loci (QTLs) 
mapping populations, using the same grandparental genotypes, with 
one from the Landes breeding population and one from a Corsican 
population, for dissecting the genetic architecture of the studied traits. 

Progeny trial 

In the 1960s, 'plus' trees were selected in the Landes de Gascogne for- 
est (southwest of France) to constitute the base breeding population 
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of maritime pine (GO generation). This study considered a half-sib 
genetic trial (HST) consisting of progenies from Gl mother trees 
selected for growth and stem straightness and crossed with a polymix 
(mix of pollen from 42 Gl trees). The HST was established in three 
sites (Escource, Lagnereau, and Cestas) with a density of 1250 trees 
ha 1 (4x2 m). A total of 196 (Escource), 156 (Lagnereau) and 174 
(Cestas) half-sib families were planted in each site, in 35 randomized 
complete blocks with a single-tree plot design. Seedlings were sown 
in 1995 and raised for 1 year in a nursery before being transferred to 
the field. These sites are characterized by a sandy podzol. Following 
Augusto et al. (2010), despite similar precipitations across the three 
sites, Escource is classified as a dry heathland, whereas Lagnereau 
and Cestas can be considered as semihumid and humid heathlands, 
respectively. These three sites correspond to those defined by Jolivet 
et al. (2007) based on understorey composition and depth of the 
water table: dry heathland (understorey vegetation dominated by 
Erica cinerea and Calluna vulgaris, deep water table), mesophylous 
(semihumid) heathland (Pteridium aquilinum, moderate depth of the 
water table), wet (humid) heathland (Molinia caerulea, water table 
close to the topsoil). There are also differences in soil nutrients, with 
less nitrogen and carbon contents at Escource compared with the 
other two sites (Supplementary Table S2 available at JXB online). 

From this progeny trial 50 half-sib families were chosen, present 
in all three sites and representing the diversity in growth observed 
at age 12 in the Gl population (Supplementary Fig. SI available at 
JXB online). Between 10 and 13 trees per half-sib family and per site 
were phenotyped for WUE (i.e. in total 634 trees at Escource, 636 
trees at Lagnereau, and 631 trees at Cestas). 

Q7L mapping populations 

A three-generation inbred mapping population (F2) was used to 
study the genetic architecture of quantitative traits. To generate the 
F2 population, selfed seeds were obtained from a single hybrid tree 
(accession HI 2) derived from a controlled cross between a Landes 
genotype (accession LI 46, from Arengosse, used as female) and a 
Corsican genotype (accession C10, from Corte, used as male). The 
closest meteorological stations (Meteo-France) were used to describe 
the climate for the origins of the parental trees: Sabres for the Landes 
genotype and Corte for the Corsican genotype. Based on data from 
1995 to 2010, Sabres showed significantly higher annual precipitation 
(1108 mm) compared with Corte (818 mm), whereas monthly mean 
potential evapotranspiration (PET) was significantly higher in Corte 
(80 mm) compared with Sabres (68 mm), especially during the growth 
period (May to August). The whole mapping population was divided 
into two seed lots. The first population (F2 greenhouse) consisted of 
200 genotypes raised in 1998 in 4-1 pots in a greenhouse and were 
subjected in their second year to a drying cycle by withholding water- 
ing for 6 weeks, after a last watering to field capacity. The second 
population (F2 plantation) consisted of 477 genotypes sown in 1998 
and raised for 1 year in a nursery before being transferred to the field. 
This trial was located at Lacanau de Mios (Gironde, France), a semi- 
humid heathland, with a density of 1250 trees ha~' (4 mx2 m) and 
trees were 9 years old at the time of phenotyping (2008). 



Phenotyping 

Traits related to WUE For the HST and the plantation populations, 
5 13 C was measured on a set of three annual rings (2006, 2007, and 
2008), which were extracted from increment cores. Extractives were 
removed by soaking the cores in pentane for 24 h. Wood was cut 
by hand into small pieces, preground in a centrifugal mill (Tecator, 
Cyclotech 1093, Sample Mill, Hoganas, Sweden) and milled to 
a fine powder in a ball mill (MM2000, Retsch, Haan, Germany). 
For the greenhouse population, 6 13 C was measured on bulk needle 
dry matter. Four pairs of needles were harvested from the upper 
third of the main stem for 200 genotypes before the drying cycle 
(5 13 C_WW_stem) and after the drying cycle (6 13 C_D_stem) for 189 
genotypes. Four pairs of needles were also sampled from a lateral 



shoot for 107 genotypes before the drying cycle (6 13 C_WW_shoot) 
and after the drying cycle (6 13 C_D_shoot) for 105 genotypes. The 
needles were ground to a fine powder using a ball mill. For all analy- 
ses, a 1 mg subsample of the plant material (wood or needle dry mat- 
ter) was combusted in an elemental analyser (Carlo Erba NA1500, 
Milano, Italy). 6 13 C was determined with a coupled continuous-flow 
isotope ratio mass spectrometer (Delta S, Finnigan MAT, Bremen, 
Germany). Carbon isotope composition (6 13 C expressed in %o) was 
calculated relative to the VPDB standard as (Craig, 1957): 

§^C= R "'~ Rsd xlOOO 

where R m and R sd are the 13 C/ 12 C ratios of the sample and the 
standard, respectively. 

For the greenhouse population, W\ was also evaluated directly 
from the ratio of net assimilation rate (A, umol irr 2 s -1 ) to stoma- 
tal conductance of water vapour (s* s , mmol nr 2 s -1 ). A and g s were 
assessed once before the drying cycle using a portable gas-exchange 
system with a cylindrical chamber (CIRAS-1 with 'pod chamber', 
PP-Systems, Herts, UK) on one pair of 1 -year-old needles per plant 
at a chamber C0 2 concentration of 1200 umol moL 1 and a satu- 
rating photosynthetic photon flux density (1500 umol m~ 2 s -1 ). The 
elevated C0 2 concentration during gas-exchange measurement does 
not affect g s in maritime pine compared with an ambient C0 2 con- 
centration (Picon et al, 1996) and was used here to allow an assess- 
ment of A under nonlimiting C0 2 concentration. The temperature 
in the measurement chamber ranged from 20 to 28 °C. Needle area 
was estimated from the length and the width of the needle included 
in the chamber, assuming a half-cylindrical needle shape. 
Traits related to growth Increment cores were sampled in the HST 
and the plantation populations. Rings were indexed by calendar year 
and ring width was evaluated using an indirect-reading X-ray densi- 
tometer (Seifert Isovolt 3003, GE Inspection Technologies, France) as 
described by Bouffier et al. (2008). In order to analyse comparable sets 
of three rings for 6 13 C and growth, this work considered ring width 
means for years 2006-2008 for the HST and plantation populations, 
corresponding to cambial age 10-12 years and 7-9 years, respectively. 

Total height and circumference at breast height were also meas- 
ured in 2008 at age 12 in the three sites of the HST. Woody biomass 
was evaluated in the greenhouse population based on weight of 
roots and primary and secondary axes without needles. 

Linkage map construction and comparative QTL mapping 

For the greenhouse population, a linkage map (here referred as 
MapF2g) was built from DNA of 200 haploid megagametophytes 
with 235 amplified fragment length polymorphisms (AFLPs), 127 
random amplified polymorphic DNA (RAPDs), and 32 protein 
markers resulting in 13 linkage groups (LGs, LG 2 was split in two 
subgroups) (Costa et al, 2000). The average distance between two 
framework loci was 15 cM (Kosambi). In this study, MapF2g was 
used in order to detect QTLs for 6 13 C, A, g s , and A/g s . 

For the plantation population, a linkage map (MapF2p) was estab- 
lished from the genotyping of 477 offspring with 248 single-nucleotide 
polymorphism (SNP) markers (C Plomion, E Chancerel, I Lesur, data 
unpublished), resulting in 12 LGs corresponding to the haploid num- 
ber of chromosomes of the pine genome. The length of the map was 
1754 cM (146 cM per LG on average) with a density of one SNP every 
7 cM. In this study, MapF2p was used to detect QTLs for 6 13 C and 
ring width. 

In order to compare QTLs detected for the greenhouse and 
plantation populations with QTLs detected by Brendel et al. 
(2002), both linkage maps (MapF2g and MapF2p) were aligned 
with two other already published linkage maps (Mapl and Map2). 
Mapl was developed in the mid 1990s with another previous set 
of haploid megagametophytes from the same three-generation 
inbred mapping population using 436 RAPD markers (Plomion 
et al., 1995) and was later completed with five single-strand 
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conformation polymorphisms (Plomion et al, 1999) and 17 pro- 
tein markers (Plomion et al, 1997). Mapl and MapF2g were 
aligned by Costa et al. (2000) based on a set of 40 common markers 
(seven proteins and 33 RAPDs). Map2 was developed by Chagne 
et al. (2002) using 620 AFLP markers genotyped on a subset of 
90 full-sib individuals (Map29 and Map2d'). The genotyping was 
extended to the whole progeny (186 full-sibs), using 219 evenly 
spaced AFLP markers segregating 1:1 (116 for the female map 
and 103 for the male map) as described by Pot et al. (2006). These 
two parental maps were used for QTL detection for 5 13 C and 
growth by Brendel et al. (2002). Mapl and Map2 were also aligned 
using 24 common markers (Single Sequence Repeats (SSRs) and 
Express Sequence Tag Polymorphisms (EST-Ps)) developed by 
Plomion et al. (1999), Chagne et al. (2003, 2004), and Guevara 
et al. (2005). Map2 and MapF2p were finally aligned by Chancerel 
et al. (2011) with 98 SNP markers. Therefore, homologous LGs 
and orientations could be unambiguously identified between the 
four maps to compare the QTLs here (MapF2g and MapF2p) 
with those of Brendel et al. (2002) (Map2). A description of the 
linkage maps is available in Supplementary Table S3 available at 
JXB online. 



Statistical analyses 

Analysis of phenotypic correlations among traits was carried out 
using linear regression analysis. 

Estimation of genetic parameters Genetic parameters were first 
estimated from the HST within each site (model 1) using the 
following linear mixed model (interactions between block and 
family effects were not significant) for univariate or bivariate 
analyses: 

y = Xb + Zf + e 

where y is the vector of phenotypic individual observations, 
b is the vector of fixed block effects, / is the vector of random 
family effects, e is the residual, and X and Z are the incidence 
matrixes linking phenotypic observations to fixed and random 
effects. 

Variances of random effects were defined as follows: 
Var(f) = G = Id.oy and Var(e) = R = Id.cr c . 

Phenotypic (op) and additive genetic (o A ) variances were obtained 
from the following equations: 

op = a} + o\ and o\ = 4 x o) 

because the families were assumed to be half-sibs. 
Genetic and phenotypic coefficients of variation (CV A and CV P ) 
were derived from these variances: 



CV A 



aAX) ,MCV P - ^ 



X 



X 



where X is the mean of the studied trait. As 5 13 C is a relative 
measurement, its coefficient of variation was only considered for 
intersite and not for intertrait comparisons (Brendel, 2014). 

Narrow-sense heritability (Ir) was estimated as: 



Ir 



The phenotypic (r P ) and additive genetic (r A ) correlations between 
two traits, y x and y 2 , were estimated as follows: 



Covp(y, y2 ) 



Cov A {y\-,yi) 



Intersite correlations were analysed from the HST, combining 
data from the three sites (Escource, Lagnereau, and Cestas) in a sim- 
ilar linear mixed model as described previously, with specific family 
covariances fitted between each site (model 2). G and R matrices 
(the variance - covariance matrix of families' effect (ie. genetic effect) 
and the variance - covariance matrix of the residual effect, respec- 
tively) can then be described as follows: 



G = 



) I t with /, being the 
number of families, and 



R 



'- °L™. J Escource © -IlMgmrem © <„,„., lam, With 

I stte being the number of individuals per site. 



Intersite genetic correlations were estimated as: 



Analyses were conducted with the restricted maximum-likelihood 
method using the software ASReml (Gilmour et al , 2009) and the 
standard errors of the estimates were calculated using ASReml 
with a standard Taylor series approximation (Gilmour et al, 2009). 
Statistical models were compared by likelihood ratio tests. 

The expected response (R) after a selection step can be expressed 
as (Lynch and Walsh, 1998): 

R = i x Op x lr 

where i-2.67 (selection of the 1% best individuals), o P is pheno- 
typic standard deviation, and h 2 is individual heritability. 
QTL mapping S 13 C, gas-exchange data (A, g s , A/g s ), ring width, and 
biomass were analysed. Data were not transformed despite slight 
deviations from normality because the interval mapping method is 
robust against this assumption (Reba'i, 1997). Moreover, permuta- 
tion analysis provides valid thresholds for nonstandard situations 
such as nonnormal distributions (Members of the Complex Trait 
Consortium, 2003). QTL detection was conducted using MultiQTL 
version 2.6 (Haifa 2005, http://www.multiqtl.com) as described by 
Marguerit et al (2012). 

As the maps of the hybrid tree were constructed using either meg- 
agametophytes (MapF2g) or needle DNA (MapF2p), QTL detec- 
tion was carried out under a back-cross model and a classical F2 
model, respectively. 

For the greenhouse population, a multi-environment analysis 
was also carried out for S 13 C estimated from the two needle posi- 
tions (stem or shoot) and sampling times before and after the dry- 
ing cycle. As no needle growth had taken place during the drought 
stress and as there were strong positive correlations between 5 13 C_ 
WW_stem, 6 13 C_D_stem, 6 13 C_WW_shoot, and 6 13 C_D_shoot, 
with slopes close to 1 (Supplementary Table S4 available at JXB 
online), these data were used as technical replicates in order to 
increase the statistical power of the QTL analysis. By reducing envi- 
ronmental variance, this method enhances the detection power of 
low-effect QTLs (Jansen et al, 1995). 



Results 

Variation in b 13 C and growth 

At 12 years (2008), trees on the three sites showed significant 
differences for growth (means, 874cm at Escource, 965 cm 
at Lagnereau, and 1043 cm at Cestas). The distribution of 
6 13 C for the HST populations (three sites) and the green- 
house and plantation populations are presented in Fig. 1 
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and Supplementary Table S5 available at JXB online. 5 13 C 
was normally distributed except for the plantation popula- 
tion (Shapiro-Wilkinson test). In the HST populations, 6 13 C 
was higher at Escource (dry site) compared to the two other 
sites, but no differences were found in CV P . As expected due 
to similar experimental conditions (field experiments with 
identical tree density, similar climate, and sandy podzol), 
6 13 C was comparable for the HST populations (mean of the 
three sites -26.1%o) and the plantation population (-26.4%o) 
but plants of the greenhouse population showed clearly lower 
6 13 C (-28.3%o) which may be attributed to the lower 6 I3 C of 
the atmospheric C0 2 in the greenhouse (mixing of outside 
C0 2 with 13 C-depleted C0 2 issued from plant and soil respi- 
ration). Both mapping populations clearly exhibited higher 
CV P for 6 13 C than the HST populations. The observed vari- 
ance was higher in the greenhouse population compared with 
the plantation population (3.1 versus 2.6%) and the green- 
house population showed a clear bimodal distribution with a 
peak ratio of 4: 1 . 

Genetic parameter estimations for traits related to 
WUE and growth 

6 13 C exhibited significant additive genetic variability in 
the three HST sites (Supplementary Table S5 available at 
JXB online). Heritability of 8 13 C was first estimated within 



each site (model 1): 0.23 at Escource, 0.41 at Lagnereau, 
and 0.34 at Cestas. These estimates (Table 1) were higher 
than heritabilities obtained for integrated growth traits (cir- 
cumference and height at age 12) but similar to heritability 
of growth (ring width 2006-2008) evaluated on the same 
annual rings as used for 6 I3 C. Because additive genetic 
variability for 5 13 C and ring width were not significantly 
different among sites, overall heritabilities of (mean±SE) 
0.29 ± 0.07 or 0.25 ± 0.07, respectively, could be estimated 
(model 2). Based on these genetic parameters, if the best 
1% of individuals were selected from the HST for 6 13 C, the 
expected response would be a shift of 0.4%o for the 6 13 C 
population mean. 

The intersite genetic correlations estimated for 5 13 C 
between Escource, Lagnereau, and Cestas with model 2 
were high (Table 2): from 0.63 between Escource and Cestas 
to 0.91 between Lagnereau and Cestas. The correlations 
between 6 13 C and growth-related traits estimated within 
each site (Table 3) were either not significant or positive. No 
significant genetic correlation was found between 6 13 C and 
circumference at age 12 whatever the site considered, but 
genetic correlations were significant between 5 13 C and height 
at age 12 at Escource (0.30, P<0.05) and Lagnereau (0.32, 
P<0.05). A highly significant positive correlation between 
6 I3 C and ring width (0.40, P<0.01) was also obtained at 
Escource. 



200 



150 



HST_Escource 
— -HST_Lagnereau 



/ / 
/ / 
/ / 

K' 



01 
01 



,_ 100 
Ol 
-Q 

E 

3 



50 



:l I 
fll 
:H 
: I I 
//' 

:/ / 
// / 
:/ / 
7 / 

/ / 



/ 



/ 



1 
1 

' \ 
\ 

I 



■HST_Cestas 
F2-greenhouse 
■ F2-plantation 



\ 



\ 



X 



-30.3 -29.9 -29.5 -29.1 -28.7 -28.3 -27.9 



-27.5 -27.1 -26.7 
6 13 C (% 0 ) 



-26.3 -25.9 -25.5 -25.1 -24.7 -24.3 -23.9 



Fig. 1. Distribution of 6 13 C in the HST populations (Escource, Lagnereau, Cestas) and the F2 greenhouse and plantation populations. 
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QTL mapping for traits related to WUE and growth 

In addition to the large variation observed for 6 13 C in the green- 
house population, a large variation for W { was also observed 
(CV P =31%). This variation was mainly determined by the varia- 
tion of stomatal conductance (g s ), which was shown by its equally 
high CV P of 31%, compared to the lower CV P for assimilation 
rate (A) of 14%. A and g s were significantly correlated (#=0.42, 
#<0.01), with a low correlation coefficient and a low slope 
(0.06 umol mmor 1 ; Supplementary Table S6 and Supplementary 
Fig. S2 available at JXB online). Further, the correlation between 
W\ and 6 13 C was significant, with a slope of 0.03, indicating a 
variation of 0.3 umol mmor 1 per %o change. Correlation between 
5 13 C_WW_shoot and A was negative and lower (#=-0.19, 
#<0.05) than that between 6 13 C_WW_shoot and g s (#=-0.53, 
P<0.01). The correlation between biomass and g s was higher 
(i?=-0.27, #<0.01) than with A (#=-0.07, not significant). 

For the greenhouse population, the multi-environment 
QTL analysis identified nine QTLs for 6 13 C (Table 4). The 
detected single-environment QTLs explained from 31% to 
59% of the phenotypic explained variance (PEV), depend- 
ing on environment (needle position and water status). The 
strongest QTL was located on LG 12, explaining 6-10% of 

Table 1 . Heritability for 6 ,3 C, circumference and height at age 
12 years, and mean ring width from 2006 to 2008 estimated in the 
half-sib genetic trial populations 

Values are mean±standard error. 





Escource 


Lagnereau 


Cestas 


6 13 C 


0.23 + 0.11 


0.41 +0.14 


0.34 + 0.13 


Circumference at 1 2 years 


0.06 + 0.08 


0.17 + 0.10 


0.27 + 0.12 


Height at 12 years 


0.07 + 0.08 


0.32 + 0.12 


0.08 + 0.09 


Ring width 2006-2008 


0.31 +0.12 


0.27 + 0.12 


0.42 + 0.14 



Table 2. Across-site correlations for b 13 C and ring width 2006- 
2008 estimated in the half-sib genetic trial populations 

Genetic correlations for 6 13 C are above the diagonal; ring width 
2006-2008 are below the diagonal. All correlations are significantly 
different from 0. 





Escource 


Lagnereau 


Cestas 


Escource 




0.89 


0.63 


Lagnereau 


0.44 




0.91 


Cestas 


0.76 


0.70 





PEV, with a QTL effect average of 0.5%o. This QTL was the 
only one detected when using a single-environment QTL 
model for each of the 6 I3 C measurements (data not shown). 
One QTL was also detected for the Alg & ratio {W-i) and g s , 
whereas none was identified for A (Supplementary Table S7 
available at JXB online). The QTLs for A/g s and g s were both 
located in the same region on LG 12, close to the QTL for 
5 13 C (Fig. 2) and explained more than 6% of PEV. 

For the plantation population, the QTL analysis identified 
one large-effect QTL on LG 12, explaining 66.7% of PEV, as 
well as a minor QTL on LG 8. The QTL on LG 12 presented 
positive additive and dominant effects (Table 5). Its confi- 
dence interval overlapped with those of the QTLs detected for 
6 13 C, A/g s , and g s in the greenhouse population (Table 4 and 
Supplementary Table S7 available at JXB online). Working 
with a genetic map of the F2 population using codominant 
markers allowed to study the contribution of alleles at the two 
flanking SNP markers of this QTL. Homozygous genotypes 
for the Landes allele as well as the heterozygous genotypes 
presented a significant lower 6 13 C compared with the homozy- 
gous genotypes for the Corsican allele (Supplementary Table 
S8 available at JXB online). The differences between geno- 
types with one or two copies of the Landes allele versus the 
homozygous genotypes for the Corsican allele (i.e. no dose of 
Landes allele) was around l%o. The interaction of these alleles 
showed that the Corsican allele was recessive. 

The phenotypic correlation between 6 13 C and ring width 
mean was not significant for the complete plantation popu- 
lation (#=-0.04) as well as when F2 offsprings of the three 
genotypes (homozygous Landes and Corsican, heterozygous) 
were analysed separately (-0.01<R<-0.05). 

Two QTLs were detected for ring width mean, explaining 
together 11% of PEV. None of them was located within the 
confidence interval of the 6 13 C QTL (Table 5). 



Discussion 

Limitations 

Quantitative genetic studies in trees are subject to the follow- 
ing limitations, which could explain the reduced heritability 
and the detection of only a few QTLs in this study. 

With respect to genetic parameter estimation, the large size 
of the HST (>5 ha for each site) implies a large uncontrolled 
heterogeneity of the soil. This can have an effect on the phe- 
notype and thus complicate the estimation of genetic param- 
eters. However, the large families used in this study made it 



Table 3. Phenotypic and genetic correlations between 5 ,3 C and growth-related traits in the half-sib genetic trial populations 



Phenotypic correlation Genetic correlation 

Escource Lagnereau Cestas Escource Lagnereau Cestas 

Circumference at 1 2 years 0.31* 0.15 0.04 0.26 0.21 -0.02 

Height at 12 years 0.31* 0.27* 0.08 0.30* 0.32* 0.08 

Ring width 2006-2008 0.44** 0.18 0.31* 0.40** 0.19 0.29* 



: P<0.05; "P<0.01 . 



The genetics of water-use efficiency and its relation to growth in maritime pine | 4763 

Table 4. Significant quantitative trait loci detected by the multienvironment analysis for water-use efficiency estimated by 6' 3 C for the F2 
greenhouse population (200 plantlets) 

Because the linkage map of the hybrid tree of the F2 mapping population was generated using haploid mega-gametophyte DNA extracted from 
each F2 plant (i.e. maternal contribution only), the QTL effect corresponds to an allelic substitution (allele a versus A). LG, linkage group; P G , 
level of significance at the genomic level; L, position of QTL on the LG estimated from original data-set; L BS , position of QTL as calculated from 
bootstrap analysis (mean±standard deviation); WW, well-watered conditions; D, drought. Asterisks indicate significant differences from zero. 
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49.6 + 26.6 


0-101.7 


0.6 
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Fig. 2. Linkage group 1 2 (LG1 2) for the four genetic linkage maps of maritime pine and associated QTLs. MapF2g corresponds to the genetic map 
published by Costa et al. (2000) from the F2 greenhouse population; Map1 corresponds to the genetic map published by Plomion ef al. (1 995) and used 
as bridge to align MapF2g to Map2o"; Map2o" corresponds to the male map published by Chagne ef al. (2002) on which SNP markers were added as 
accessory markers by Chancerel ef al. (201 1). The framework Map2 was used for QTL detection of water-use efficiency by Brendel ef al. (2002); and 
MapF2p corresponds to the map established by Chancerel ef al. (201 3) from the analysis of the F2 plantation population. Bars indicate the range of 
mean positions from L and L BS (see Tables 4, 5 and Supplemental Table S7 available at JXB online). Whiskers are the 95% confidence interval around L BS 
(this figure is available in colour at JXB online). 



possible to establish trials consisting of blocks, which can be 
used to take this heterogeneity into account. 

Many tree species cannot be easily propagated vegetatively 
(including maritime pine), therefore QTL studies in these spe- 
cies are based on a single-tree plot design without 'clonal' rep- 
licates. Consequently, environmental heterogeneity, especially 
of the soil, cannot be controlled by introducing blocks into the 
experimental design. Spatial variation of soil properties could 
therefore have an impact on the estimation of the phenotype 
and thus result in less statistical power to detect low-effect 



QTLs. Therefore, before the data analysis, this work verified 
for the QTL field experiment, using a spatial analysis of the 
phenotype of the measured trees, that no spatial pattern was 
detectable (Supplementary Fig. S3 available at JXB online). 

Traits of interest (growth, wood properties, resistance to 
biotic and abiotic stresses) often represent an integration over 
several growing seasons, which will average out annual envi- 
ronmental variation when the genotype x environment inter- 
action is not particularly strong and thus result in a more 
robust estimation of phenotypic values for each genotype. 
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Table 5. Significant quantitative trait loci of water-use efficiency estimated by b 13 C and ring width mean for the F2 plantation population 
(trees) 

LG, linkage group; P G , level of significance at the genomic level; L, position of QTL on the LG estimated from original data-set; L BS , position of 
QTL as calculated from bootstrap analysis (mean±standard deviation); PEV, percentage of explained phenotypic variance; PEV a , percentage of 
explained additive variance; a, additive allelic effect; d, dominance allelic effect. Asterisks indicate significant differences from zero. 
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a 


d 


5 13 C 


447 


8 


4.4 


0.013 


60.0 


47.9 + 20.7 


7.2-88.5 


1.6 


0.1 


-0.04 


-0.18 






12 


81.5 


0.0005 


6.0 


6.1 +0.9 


4.4-7.8 


66.7 


46.3 


1.38* 


-0.65* 


Ring width mean 


460 


1 


5.2 


0.0008 


142.2 


137.1 +20.9 


96.1-144.8 


4.9 


3.3 


-0.81* 


0.39* 






10 


6.9 


0.0006 


12.9 


23.7 + 24.3 


0-71 .3 


6.1 


5.6 


1.05* 


0.22 



Genetic parameters of the breeding population 

Genetic parameters were estimated in the Landes provenance 
as this genetic material constitutes the main French mari- 
time pine breeding population. 5 13 C heritability measured 
across different environments was moderate (0.29 ±0.07) 
and did not vary significantly among sites. It was similar 
to that found for ring width growth (0.25 ±0.07) evaluated 
for the same years but slightly higher than for total height 
and circumference at the age of 12. The genetic parameters 
of 5 13 C were accurately estimated in this study due to the 
presence of a large genetic variability in the Gl breeding 
population on the one hand and the analysis of the same 
genetic material in three sites contrasting in water availability 
on the other. Thus, this robust estimation of 6 13 C heritabil- 
ity was obtained using a wider genetic background within 
the breeding population compared to the first estimation of 
0.17 ±0.06 provided by Brendel et al. (2002) using a 12 x 12 
half-diallel (i.e. based on 12 genotypes). In the current data 
set of more than 1800 trees, the genetic variability for 5 13 C 
was as large as the one found across the natural range of the 
species (Guehl et al., 1996). Despite the different soil con- 
ditions, the genetic within-site variability was not different 
among sites, suggesting only a few genotypexenvironment 
interactions. 

In other coniferous species, the highest narrow-sense herit- 
abilities estimated for 6 13 C were reported by Prasolova et al. 
(2001) for Araucaria cunninghamii (0.66-0.72) and by Johnsen 
et al. (1999) for Picea mariana (0.54), whereas very low values 
(<0.1) were observed for Pinus taeda (Baltunis et al., 2008). 
These contrasting estimates illustrate the complex determin- 
ism of 6 13 C and its underlying leaf-level processes which are 
under both genetic and environmental control. On a more 
technical side, this variability in heritability estimates can also 
result from poor experimental designs in terms of sample size: 
22 and 26 open-pollinated families with 101 and 130 trees for 
Prasolova et al. (2001); 18 full-sib families (from seven parents 
only) with 820 trees for Johnsen et al. (1999); 61 full-sib fami- 
lies (from 32 parents) with 1027 trees for Baltunis et al. (2008). 

Intersite correlations for 5 13 C revealed significant but 
low genotypexenvironment interactions, which suggest that 
the ranking of families was conserved among sites. Sites 
that were environmentally closer showed a stronger corre- 
lation, suggesting that genotype-environment interactions 
might become stronger when more contrasting planting 



environments are considered. Moreover, mostly the inter- 
site correlations for 6 13 C were higher than those obtained 
for growth traits, indicating stronger genotypexenvironment 
interactions for the latter. This finding is interesting from a 
breeding perspective, as selection for higher WUE might be 
more stable than selection for growth. 

Genetic architecture of 8 13 C and W, in a Landes x 
Corsican genetic background 

The genetic architecture of 6 13 C has been previously studied 
within the Landes provenance genetic background (Brendel 
et al. , 2002) and one strong QTL for 6 13 C was detected. Here, 
a wider genetic background, a Landes x Corsica F2 mapping 
pedigree, was used for QTL detection for WUE- and growth- 
related traits in two independent experimental settings: a 
greenhouse with 2-year-old seedlings and a field plantation 
with 9-year-old trees. 

In the plantation population, a major QTL region linked 
to 6 I3 C on LG 12 was detected (accounting for 60% of PEV) 
that was also detected in the greenhouse population. This 
QTL had not been detected in the intra-Landes pedigree 
used by Brendel et al. (2002), and vice versa, the major QTL 
detected by Brendel et al. (2002) was not detected in this 
study (Supplementary Table S9 available at JXB online). 
This suggests that the underlying genetic factors are dif- 
ferent between the pedigrees. Compared to the geographi- 
cal origin of the Landes parent, the original population of 
the Corsican parent is subjected to higher annual potential 
evapotranspiration, especially in the summer. This results 
probably in a higher summer water deficit, which might have 
caused specific strategies (therefore different genetic mecha- 
nisms) to cope with soil water deficit. 

The greenhouse experiment allowed a more detailed study 
of traits related to WUE. Here, the alignment of the differ- 
ent LGs (Fig. 2) allowed this work to show that the major 
QTL for 6 13 C located on LG 12 in the plantation popula- 
tion: (i) colocalized with the QTL detected in the greenhouse 
population; and (ii) colocalized with a QTL for a direct 
measurement of W\ and with a QTL for g s . Similarly to the 
Brendel et al. (2008) study on Quercus robur L., the colocali- 
zation of the QTL as well as the positive phenotypic correla- 
tion between 6 13 C and A/g s confirmed the major QTL for 
5 13 C as a QTL for WUE rather than for one of the Farquhar 
et al., (1989) model parameters, which depend on biology 
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(e.g. day respiration, photorespiration, mesophyll conduct- 
ance of C0 2 ). Further, the colocalization with a QTL for g s 
suggests an introduction of alleles by the Corsican, suppos- 
edly more drought adapted parent, which might be related 
to stomatal responses, resulting in genotypes that are more 
efficient in water use. The correlations between gas exchange 
related traits and 6 I3 C suggest that the variability of WUE 
in this cross was probably driven by variation in g s rather 
than A. The coefficients of variations of A and g s suggest a 
larger diversity of the latter in the greenhouse population. 
Additionally, the significant correlation between the two 
traits had a low slope, suggesting only a relatively low impact 
of g s on A. Thus, within this family, a lower g s effectively 
reduced water consumption, increasing WUE with only a 
low impact on aboveground biomass growth. It is worth not- 
ing that the impact of this genomic region on LG 12, which 
controlled the variability of 5 13 C, became much stronger in 
the field plantation (8-67%), even though these trees were 
much older than the seedlings in the greenhouse, and meas- 
urements were done on wood samples in the field, which inte- 
grated whole-canopy WUE as well as post-photosynthetic 
discrimination (Badeck et al, 2005). This suggests that the 
differences in stomatal conductance among the seedlings 
were probably cumulative over time and the crown, resulting 
in a much stronger differentiation among older trees and no 
detectable relationship with aboveground growth. As far as is 
known, this is the first time that such a strong QTL for WUE 
estimated by 8 13 C in a perennial plant has been identified. 
Similar strong-effect QTLs (21-64% of PEV) have only been 
detected in recombinant inbred lines on Arabidopsis thaliana 
(Masle et al, 2005) and Solarium spp. introgression lines (Xu 
et al, 2008). For a perennial species, Q. robur, a rather strong 
QTL (21-31% of PEV) for 6 13 C has also been detected with 
a multiyear QTL analysis (Brendel et al, 2008). In other 
perennial species, 2-17 QTLs for 8 13 C have been detected, 
explaining each from 2 to 18% of PEV: Pinus pinaster 
(Brendel et al, 2002), Gossypium interspecific cross (Saranga 
et al, 2004), Castanea sativa Mill (Casasoli et al, 2004), 
Salix interspecific population (Ronnberg-Wastljung et al, 
2005; Weih et al, 2006), Populus interspecific cross (Street 
et al, 2006; Monclus et al, 2012), and woody fruit crops 
( Vitis spp, Marguerit et al , 2012; Mains domestica, Regnard 
et al, 2009; Segura et al, 2009). In these studies, sample 
sizes varied from 92 to 331, which, according to simulations 
by Beavis et al. (1994), would suggest that QTL effects were 
overestimated. In annual crops, genetic control of 8 I3 C has 
been shown to be under the control of several independent 
loci explaining a small part of PEV, not repeatable across 
years, and thus subject to large QTL-environment interac- 
tions (Rebetzke et al. 2008). However, this is not necessarily 
the case, as shown by Brendel et al (2008) for pedonculate 
oak, where several QTLs were detectable across years. In the 
present study, the relevant choice of parents with contrasting 
behaviour under water deficit (Landes and Corsican geno- 
types), the large variability present in the F2 population, a 
potential cumulative effect with ageing for the QTLs, and the 
size of the population (447 phenotyped trees) allowed the 
reliable detection of a large-effect QTL. 



Relationship between 8 C and growth traits 

At the interprovenance level, several studies have analysed 
the correlation between 6 13 C and growth in maritime pine. 
Nguyen-Queyrens et al (1998) found a positive phenotypic 
correlation between tree height and 6 13 C with trees from the 
French Landes and the Moroccan Tamjoute provenance. On 
the contrary, Correia et al. (2008) found a negative correla- 
tion between height and 6 13 C across provenances from the 
Atlantic coast (France, Spain, and Portugal), as did Corcuera 
et al. (2010) across three Spanish populations and Lamy et al 
(201 1 ) for French, Spanish, and Moroccan provenances. 

At the within-provenance level, Brendel et al. (2002) 
showed a positive phenotypic correlation between ring width 
and 6 13 C but this correlation appears not significant at the 
genetic level. Enlarging the literature review to other for- 
est trees evaluated in field after the juvenile age, no general 
trend could be highlighted. Some authors found negative cor- 
relations, such as Xu et al. (2000) in Pinus elliottii x Pinus 
caribaea, Johnsen et al. (1999) in Picea mariana, and Lauteri 
et al. (2004) in C. sativa, while others showed positive cor- 
relations, such as Prasolova et al. (2001) and Xu et al (2003) 
in Araucaria cunninghamii or nonsignificant correlations as 
Ponton et al. (2001) in Quercus petraea. These results can be 
first explained by a poor estimation of the genetic correla- 
tion due to limited sampling, as 6 13 C measurements are cost 
and time demanding. This discrepancy also suggests that the 
direction of this correlation is driven by different physiologi- 
cal mechanisms depending on the species and even ecotype 
considered. A positive correlation between WUE and growth 
could suggest that WUE is primarily determined by leaf 
assimilation, which would impact variability of biomass 
accumulation. On the other hand, a negative correlation 
could indicate that WUE is primarily determined by stomatal 
conductance, having little impact on biomass accumulation. 

No significant unfavourable genetic correlations were 
found between 8 13 C and growth traits. Thus, it appears possi- 
ble to select maritime pine for 8 13 C without counter selection 
for growth within the Landes genepool. Moreover, highly 
significant positive (favourable) correlations were estimated 
from the dry site (Escource), suggesting that the acclimation 
to this site resulted in an increased WUE that was genetically 
linked with an increase of growth. 

As all three planting sites are within a small region in the 
southwest of France, the overall climate is similar and the 
differences in soil water content are mainly dependent on the 
soil type. The driest soil at Escource, which contained the 
highest proportion of sand, also showed relatively low nitro- 
gen and phosphorous contents. The resulting nutrient limi- 
tation was therefore probably higher in the dry site due to 
reduced nitrogen mineralization compared to the other sites 
(Stanford and Epstein, 1974; Cassman and Munns, 1980). 
Thus, in addition to variability in WUE among genotypes, 
variability in nutrient-use efficiency, which has an impact on 
photosynthetic capacity, could be the cause of the observed 
positive genetic correlation between growth and 6 13 C. 
Disentangling the effect of both factors will require further 
investigation. 
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In the plantation population, there was no coincidence 
between the major QTLs for 5 13 C and growth traits. In addi- 
tion, there was no phenotypic correlation (R--0.04, not sig- 
nificant) between these two traits. Additionally, the test of the 
homogeneous group at the closest markers from the QTLs did 
not show the same pattern of significant differences between the 
homozygous and heterozygous groups (Supplementary Table S8 
available at JXB online), suggesting no common genetic regula- 
tion between growth and S 13 C. These results therefore suggest 
that the genetic determinisms of 5 13 C and ring width are inde- 
pendent. The absence of relationships between growth and 6 13 C 
has also been found in Populus interspecific families (Marron 
et al. , 2005; Monclus et al, 2005, 2006, 2012; Dillen et al. , 2008). 
The absence of a trade-off between these traits suggests that it 
should be possible to select maritime pine genotypes combining 
high growth and high WUE. 

Conclusions 

This work presents a comprehensive genetic study over multi- 
ple sites and different ages that quantify the genetic relationship 
between growth and WUE. A large sampling of the maritime 
pine breeding population (Landes provenance), established in 
sites contrasting in water availability, showed medium heritability 
for 6 13 C, an indirect measure of intrinsic WUE. No unfavourable 
correlation with growth was found in this population. The genetic 
variability for 6 13 C appeared higher in a Landes x Corsica F2 
population and the genetic architecture revealed a strong QTL 
detected in the same genomic region in two experimental settings. 

The expected genetic gain for a selection for 5 13 C in the 
Landes breeding population was rather low (shift of 0.4%o if 
the best 1% of genotypes are selected). However, the Corsican 
genetic background, which brings a recessive allele improving 
WUE greatly by decreasing stomatal conductance with only 
little effect on assimilation rate and growth, could be used in 
controlled crosses with the Landes genetic background. More 
generally, it seems interesting to explore the genetic diversity 
for WUE within the Corsican gene pool for implementation 
in the maritime pine breeding programme. The usefulness of 
WUE, estimated by 6 13 C, is under debate in breeding of annual 
crops for higher yield under drought conditions (Blum, 2009; 
Chen et al., 201 1; Sinclair, 2012), with the main argument that 
often higher WUE is linked to lower stomatal conductance, 
which then also decreases carbon assimilation and thus yield; 
however, for the maritime pine genetic backgrounds studied 
here, higher WUE was not genetically linked to less growth 
and there is therefore space for improvement in water use. 

The genetic improvement of the Landes breeding popu- 
lation in itself for WUE might not be a major aim. Rather, 
more research is needed to support the interest of the Corsican 
genetic background in hybrid combination with the Landes 
provenance, which is well adapted to the Aquitaine region, 
especially in terms of frost resistance. To gain more insights 
into the functional background of observed variability in 
WUE, an approach of choosing specific genotypes for detailed 
studies, as has been done for pedonculate oak (Brendel et al. , 
2008; Roussel et al, 2009a, b), might be feasible. The genetic 
material created within the maritime pine breeding programme 



(recurrent selection scheme between the Landes and Corsican 
genepools in the early 1980s and large progeny provenance tests 
of Corsican ecotypes in the 2010s) provide further research 
grounds to study the genetics of WUE in this species. 
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